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In preparation for an experimental study of magnetorotational instability (MRI) in liquid metal, 
we explore Couette flows having height comparable to the gap between cylinders, centrifugally stable 
rotation, and high Reynolds number. Experiments in water are compared with numerical simula- 
tions. Simulations show that endcaps corotating with the outer cylinder drive a strong poloidal 
circulation that redistributes angular momentum. Predicted azimuthal flow profiles agree well with 
experimental measurements. Spin-down times scale with Reynolds number as expected for laminar 
Ekman circulation; extrapolation from two-dimensional simulations at Re < 3200 agrees remarkably 
well with experiment at Re ~ 10 6 . This suggests that turbulence does not dominate the effective 
viscosity. Further detailed numerical studies reveal a strong radially inward flow near both endcaps. 
After turning vertically along the inner cylinder, these flows converge at the midplane and depart 
the boundary in a radial jet. To minimize this circulation in the MRI experiment, endcaps con- 
sisting of multiple, differentially rotating rings are proposed. Simulations predict that an adequate 
approximation to the ideal Couette profile can be obtained with a few rings. 



I. INTRODUCTION 

Laboratory experiments using liquid metal have made important tests of magnetohydrodynamic (MHD) theories and 
processes. Examples include laboratory demonstration of the Alfven wavep], magnetofluid convection 0, and dynamo 
actionQ- Recently, a new kind of liquid-metal MHD experiment, motivated by astrophysics, has been proposed by 
the authors0,|5| to study magnetorotational instability (MRI), which is believed to dominate the transport of angular 
momentum in electrically conducting accretion disks. Liquid gallium will be used in a Couette flow between cylinders 
of radii n < r2 and angular velocities < fl2 < f2i but > rffli so that the flow will be stable against conventional 
Taylor-Couette instabilities (TCI). 

Both TCI and MRI are governed by the radial profile of azimuthal velocity, v v (r). In an inviscid fluid, TCI occurs 
wherever the specific angular momentum, J = rv v , decreases outwards, dJ 2 /dr < 0. If the fluid is also a perfect 
electrical conductor, MRI occurs when the angular velocity, £1 = v v /r, decreases outwards, <Kl 2 jdr < 0. In this case, 
MRI can occur in an arbitrarily weak axial magnetic field; the field must at any rate be weak enough so that the 
transit time of Alfven waves across the flow is less than f2 _1 . The usual situation in accretion disks, which are often 
excellent conductors, is fl cx r~ 3 / 2 (Kepler's law) so that TCI is stabilized but MRI is not. In a fluid with large 
resistivity and finite but small viscosity, such as liquid metals, the range of flow parameters unstable to MRI shrinks 
significantly, while the range unstable to TCI is effectively unchanged. Therefore, laboratory flows must be set up 
with precision in order to demonstrate MRI while suppressing TCI. 
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In previous linear stability analyses of gallium Couette flow|j,|5j], we adopted periodic boundary condition in the 
vertical (axial) direction, ignoring the effects of the top and bottom interior surfaces of the vessel ("endcaps"). The 
choice of vertical boundary conditions is probably inconsequential when the height of the flow (H) is much larger 
than the gap width, as in Taylor's classic experiments [||. Our experimental volume 7ri?(r| — rf) will be limited by 
the availability of gallium, a far more expensive fluid than water, while the gap must be wide enough so that the 
magnetic diffusion time is not much shorter than the rotation period. These considerations drive us to an aspect ratio 
H/(r2 — ri) ~ O(l), in which the endcaps may assume great importance. 

We have performed a water experiment and complementary numerical simulations to study the effects of the endcaps 
and, if possible, to find a way to set up a short Couette flow that is unstable to MRI yet stable against TCI. Since the 
viscosities of the two fluids are similar, standard visualization techniques in water serve to predict the flow structure 
in the opaque liquid gallium, at least in the absence of magnetic field. 

Since the pionee ring work b y Benja min[ fl |a>l p| , TCI in finite size cylinders have been studied in detail with stationary 
outer cylinders El El El El M, El E OMllllIi, and with rotating outer cylinders [2l|. In our TCI-stable flows, 
the outer cylinder must rotate. To allow a wide gap, we use a relatively small radial aspect ratio 77 = r\/n = 0.256, 
so that the commonly used narrow-gap approximation does not apply. The rotation rate of the inner cylinder (Oi) is 
so high that the Reynolds number 

Jte^ ri(ra ; ri)ni ~Q(10 6 ), (I) 

is orders of magnitude larger than typical Couette flow experiments. In connection with research on the tropospheric 
jet stream, Dunst performed a water experiment in a short cylindrical annulus with parameters similar to ours|22|. 
We will comment on this paper in section Ivl 

The outline of the paper is as follows. Sec|H] describes the experimental apparatus and most of our experimental 
results, except for spin-down measurements, which are deferred to Sec II VI Numerical methods are described in Sec lIIII 
and numerical simulations are presented and compared with experiment in Scc lIVI SecE] contains a summary of our 
main results and a discussion of their significance. 

II. WATER EXPERIMENT 

A. Experimental Apparatus 

The experiments are performed in a short circular Couette flow illustrated in Fig. ^ A container, made of trans- 
parent acrylic plates and cylinders, is mounted on a stainless steel flange which is driven by a DC motor (1.25 HP, 
by Creative Industries). The speed is measured by a laser-based tachometer. The inner cylinder, made of aluminum, 
is inserted through a lip seal from the top and is driven in the direction of the container by an AC motor (3 HP, 
by Lincoln Motors) with a controller. For given speeds of inner and outer cylinders, f2i and f^, the flow typically 
requires about 10-30 seconds to reach a steady state, in which most of the flow measurements are performed. 

B. Measurement of Azimuthal Flow Profiles 

In order to measure flow profiles as a function of radius and height, small particles with sizes on the order of 1 /im 
made of mica and titanium dioxide ("Sparkle"by Lee S. McDonald, Inc.) are mixed into the working fluid (water). 
Being small, the particles follow the flow well. A sheet of light of approximate thickness of 0.5 cm is generated by 
a horizontal slot in front of a bright halogen light source. An electronically gated intensified-CCD camera (by ITT 
Corp.) images the flow illuminated by the sheet, whose height is adjustable. Images are saved to a PC using a frame- 
grabber at a rate of 60 images per second. Particles appear in the images as streaks, whose length indicates flow 
speed. By combining measurements at many radii and heights, the azimuthal flow can be mapped out as a function 
of r and z. The measurements were calibrated, at various heights, by imposing uniform rotation, f2i = O2 = 150 rpm. 

The experimental results are shown in Fig. [3] for the case of f^i = 2000 rpm and ^2= 150 rpm. There are a few 
characteristics worth mentioning here. First, the measured velocity is significantly smaller than in an ideal, infinitely 
long Couette flow having the same £li and [ca. ljAlfl ]. The discrepancy is largest at small radii. Secondly, the 
velocity must rise sharply from ~ 3 m/s to match the inner cylinder at ~ 8 m/s. Unfortunately, diagnostic access to 
the flow is limited near the inner, outer, and top boundaries. Thirdly, the azimuthal velocity decreases with radius 
except at the locations near the outer edge, whereas it would decrease everywhere in the ideal Couette flow. Fourthly, 
the dependence on z is at most comparable to the experimental errors. This is consistent with Taylor-Proudman 
theorem |23l | which predicts small z variations in a rotating flow with small viscosity. 



3 



ICCD 



n 



Mgtor 



Acrylic 



Light Source 



Light Sheet 



Aluminum 
Rl 



(tot or 



FIG. 1: Schematic of experimental apparatus. Radii of the inner and outer cylinders are ri=3.8 cm and r2 = 14.9 cm, 
respectively. The height is H=1Q cm. 
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FIG. 2: Measured v v profile at different z when Qi — 2000 rpm and ^2= 150 rpm. 



The observed profile of azimuthal flow has unfavorable implications for the proposed MRI experiments. The goal 
is to set up a flow unstable to MRI while stable to TCI. However, the sharp decrease of v v near the inner cylinder 
will certainly incite TCI while the rest of the flow, because v v falls more slowly than intended, will be more resistant 
to the MRI. As a result, the system as a whole could have a mixture of both instabilities or, even worse, only the 
TCI. The observed deviations from ideal Couette flow are due to the endcaps. They need to be understood and, if 
possible, to be minimized in order to demonstrate MRI unambiguously. We note that a proposed MRI experiment 
using sodium [24| should suffer from the same complications since its aspect ratio is also small, viz. H/(r2 — ri) = 2. 
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III. NUMERICAL METHODS 

A. Mathematical Model and Algorithm 

We have developed a 2-dimensional computer code to simulate the water experiment described in the previous 
section. We use the stream function-vorticity method[2sfl in cylindrical coordinates (r, ip, z), as illustrated in Fig. 
Assuming axisymmetry, <9 v (tv, v Vl v z ) = 0, and incompressibility, V • v = 0, we can express the poloidal components 
of the flow in terms of a stream function, ip: 

1 di> 

Vr = (2) 
r oz 

1 dip 

v z = — -5-. 3 
r or 



Z^-—— + ZT^2)^ = UJ v ( 4 ) 



The azimuthal component of vorticity becomes 

'Id 2 Id Id 2 
s r dr 2 r 2 dr r dz 2 

Since the normal component of velocity vanishes, the boundary conditions for ip are 

ip = at r = ri,r2, (5) 
ip = at z = ±H/2. (6) 

The boundary conditions for v v are as follows. At the inner and outer cylindrical walls, 

v v = rxOx at r = r%. (7) 
v v = r 2 VL 2 at r = r 2 . (8) 

Since the top and the bottom endcaps of the vessel rotate with the outer cylinder in our apparatus, the boundary 
condition there are 

v v = —fl 2 at z = ±H/2. (9) 
ri 

The no-slip conditions on v r at z = ±H/2 and on v z at r = r\,r% yield boundary conditions for u v via eqs. (^Ql : 

1 d 2 ip 

= r~gZ2 at r = r i: r 2 I (10) 
1 d 2 ib 

= -ITT at z = ±H/2. (11) 

r az z 
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The fundamental variables in our numerical simulation are v v and lo v . Their governing equations are 

du) w dD r dD, /_„ 1 . 

~W = -af-^ + v { v -^H' (12) 

-w = ~V^ + v ^ + —) +v V (13) 



where 



J /" 9J^\ i oj v ( dij) 



r r 2 \ dr ) r \ <9r ' 
uo v ( d%jj\ t J ( dJ 



r 



dz I r 2 \ d 



D z = ^ +^ : (15) 



J = rw^, (16) 



and 



<9 2 1 9 d 2 
5r 2 r <9r 9z 2 



(17) 



We use second-order spatial differences on a uniform grid, with typical size N r xN z = 100 x 100, and a fourth-order 
Runge-Kutta method for the temporal integration, with typical time step St = 7.273 x 10~ 3 (r2 — J"i)/rifii. 
The algorithm for each time step is as follows. 

1. Integrate the basic eas. (ll2ll and (jlHfl to get new aj v (i,k) and v v (i,k) in the bulk region (2 < i < N r — 1, 
2 < k < N z - 1). 

2. Solve eq.Q for ^ with the new lo v as source term and boundary conditions JSJ-JHI). 

3. Set the boundary values of uj v and v v using eqs.(|7J), (JSJ, lfTU|) and (fT*Tl> . 

4. Get the auxiliary variables u r , u z , D z , and J from eqs.©-© & (|14fl - i|16[) . 



B. Simulation Method 



The boundary conditions 0, JHJ, and © require a jump of v v (r) at the corners (r — r\, z — ±H/2). A commonly 
used technique to avoid this singularity is to make a small "buffer region" r\ < r < n + e, in which v(r,±H/2) 
varies smoothly 26] . But we found that such a "buffer region" is not necessary in our scheme. Following the standard 
technique, we first set e to 10% of the gap width r-i — r%, and then gradually reduced it to the radial mesh size, Ar. 
In the latter case the angular velocity of the boundary jumps from f2i to between the innermost two grid points 
(i = 1,2). We confirmed that simulation results are not affected by the size of e. All calculations shown in this paper 
use e = Ar. 

All calculations begin with both the boundaries and the fluid at rest. The boundaries (side walls and endcaps) 
accelerate to their final angular velocities over a short time interval < t < r. We confirmed through many trials that 
the final state of the flow does not depend upon r. All simulations shown here were calculated with r = (ri — r2)/(rif2i), 
our unit of time. 



IV. RESULTS OF SIMULATIONS 



A. Parameters and units 



The simulations shown in this section use the same dimensions and rotation rates as as the experiments described 
in Sec.[n] The computational units of length, velocity, and time are I = — ri, v = rifli, and £/v. 

The Reynolds number of the experiment [eq.Q] is too high for direct numerical simulation. So we start by 
simulating very low Re and increase it until we find characteristics of the flow that depend only weakly on Re or 
follow a clear scaling. The largest simulated Re that we report is 3200. 



FIG. 4: Simulated v v profile when Re = 1. 



B. 



Low Re flows 



The sudden acceleration of the boundaries at the beginning of the simulation induces an intense flow inside the 
vessel. At small values of the final rotation rate, i.e. small Re, the fluid quickly relaxes to a laminar steady state. 

Figure 21 shows the azimuthal speed, v v , for Re = 1. The maximum speed is achieved at the inner cylinder (on 
left in Fig. 0J, where v v = 1. In this low Re limit, poloidal flow is almost absent; the maximum value of v r and v z 
is 7.5 x 10~ 4 . Note that the jump in Vl at the lefthand corners of the boundary does not prevent a smooth solution 
elsewhere. 

Actually, as Re — > 0, the solution for the flow can be found in closed form. If we set v r — v z = lu v = in the basic 
equations (|12fl - 1|13|) . the stationary azimuthal flow v v satisfies 



with the boundary conditions 10-10 • Wendl[22| has given the analytical solution of this equation with slightly 
different boundary conditions, corresponding to Q 2 = 0. The solution to our problem is given simply by adding a 
uniformly rotating component rf2 2 to Wendl's solution. The v v profile shown in Fig. 0] is essentially identical to the 
analytical solution thus constructed. This serves as one benchmark for our code. 



We have seen that poloidal flow is negligible in the low Re regime. As we increase Re, poloidal circulation develops 
and the azimuthal flow changes. Figure [5] shows cross sections in the poloidal (r, z) plane for different Re numbers, 
from 100 to 3200. Each panel in Figgis a snapshot of the nonlinearly saturated state. 

Figure EJa) shows that v v has much the same pattern at Re — 100 as at Re — 1 (Fig. HJ. In both cases, it is 
symmetric about the horizontal plane z = and time independent. The flow becomes unsteady at Re > 400. The 
asymmetric profile of Re = 800 [Fig.OJd)] results from unsteady flow. All flows above this Reynolds number fluctuate, 
with an amplitude that increases with Re. 

One of the important features shown in Fig. [5] is that the contours of v v tend to be parallel to the rotation axis 
{£). This is a manifestation of the Taylor-Proudman theorem, viz., that low- frequency horizontal motions tend to be 
independent of height in an inviscid fluid rotating about a vertical axis|28|. The Taylor-Proudman theorem is usually 
discussed for an almost rigidly rotating fluid, but as shown in the Appendix, a similar tendency exists in differentially 
rotating flow provided dJ 2 jdr > 0. 

Another characteristic feature of the v v contours in Fig.JSJis their tendency to concentrate towards the inner cylinder 
(at left) at large Re. This shows the development of a boundary layer. Boundary layers also develop on the top and 




(18) 



C. High Re flows 
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bottom endcaps [see panels (e) and (f)]. Note also the increasingly sharp protrusion of the contours on the inner 
cylinder near the midplane. 

The development of the poloidal (v r ,v z ) flow is equally interesting. It is shown in the vector plots of Fig. for 
the same simulations as in Fig. |SJ Vectors appear at every third grid point in r and z. Their lengths indicate that 
the poloidal component becomes stronger at higher Re. The arrow lengths are linearly normalized. The amplitude 
of the unit velocity, that is the rotation speed of the inner cylinder, is indicated by arrows on the top of the panels 
(a) to (f). A radially outward, jet-like flow is seen near the midplane in the higher- Re simulations. It is undoubtedly 
the counterpart of the spike in the v v contours in Fig. [S] This jet-like flow was not expected before we began our 
numerical simulations. Its structure will be analyzed later. Here we show the temporal behavior of the poloidal flow, 
including the jet. 

The jet flaps unsteadily. Figure [7| shows the time development of the stream function, ip, for Re = 3200 over 
almost half of the period of oscillation of the jet. The period of the oscillation is about 36.95 in the normalized time. 
Figure Ufa) is a snapshot taken at t = 1326.68, and subsequent panels (b) to (f) are taken with regular 5.82 interval 
time. Throughout the oscillation, the root of the jet (at r = r\,z = 0) remains fixed while its tip flaps violently up 
and down. The vector plot in panel (f) of Fig. shows a snapshot of this flapping motion at t = 1294.67. 



Figure [S] shows v v (r) for different Reynolds numbers (200 < Re < 3200) at z — 4cm. (This is the height where the 
most reliable experimental data can be obtained.) One can see from this figure that the v v profiles are consistent for 
the higher Re flows, Re = 800, 1600, and 1600. Note also that the profiles for Re = 1600 and 3200 are very close, 
suggesting possible convergence at large Re. 

In Fig.0 we have superimposed the experimental data for v v at z = 4cm (shown also in Fig. EJ) on the corresponding 
profiles from simulations at Re = 1600 and 3200. The agreement is remarkable when one considers that the highest 
Re achieved in simulations, which are in 2D, is a factor of 300 smaller than that of the experiment, which is 3D. 

We have seen in Fig. that the high Re (> 400) flows are time dependent. To show the amplitude of the temporal 
fluctuations, snapshots of v v for Re = 1600 are superimposed in Fig. EI] (The curves in Fig. are time averages.) 
We can see that the amplitude is larger in the outer half of the flow (r > 9cm). This can be explained by the flapping 
motion of the jet (see Fig. EJ. 

Both experiment and simulation indicate that the v v profile is concave in the outer half of the flow (Fig. and 
that u„(r) is an increasing function of radius near the outer cylinder (12cm < r < 15cm). The fluid in this region 
rotates almost rigidly, on average, at the angular velocity of the outer cylinder and endcaps (^2)- 

Snapshots of the angular momentum J = rv v are shown in Fig. llll for the same simulation and as in Fig. 1101 Note 
that dJ/dr > for the ideal Couette flow (solid curve) since we aim for stability against the Taylor-Couette mode. 
An interesting feature of this figure is that the simulated J curve is even flatter than the ideal profile. This is a 
consequence of the poloidal circulation and jet, which tend to mix angular momentum in the interior of the flow. 



In order to extract essential features of the spatial structure of the flow, it would be convenient if we could suppress 
its temporal fluctuations, especially the flapping of the jet. For this purpose, we have made a set of simulations with 
the following symmetries imposed: 



Fig. ^] shows the result when Re = 1600. One can perhaps interpret this state as an average of the flow over one 
period of the jet's flapping motion. 

The panel (a) in Fig. 1121 clearly shows the structure of the poloidal circulation. It consists of three main parts; (i) 
inward flow in boundary layers at the endcaps; (ii) axial flow towards the midplane on the inner cylinder; and (iii) an 
outward jet centered on the midplane. 

The boundary layer on the inner cylinder has a characteristic triangular shape [panel (b) in Fig. I12j . The contours 
are "squeezed" by the inward boundary-layer flows near the top and the bottom caps, while, the tip of the triangle is 
"pulled" by the outward jet flow. 



D. Profile of Azimuthal Flow 



E. Jet 
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FIG. 5: V v profile. 






FIG. 6: Poloidal flow profile. 
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FIG. 7: Contour plots of the stream function in a time sequence showing the oscillatory motion of jet (Re = 3200). 
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FIG. 8: Simulated v v profile at z = 0. 



The contour lines of v v and SI [Fig. I12f and (c)] suggest that the flow attempts to obey the Taylor-Proudman 
theorem outside of the boundary layers and the jet. Since the latter regions are thin, especially at high Re, viscous 
forces becomes important there and the Taylor-Proudman theorem does not hold. 

The contour lines of O, Fig. I12f c'l. show that the the outer half of the fluid rotates approximately at the angular 
velocity of the outer cylinder and endcaps, as we have already seen in Figs. HOIandllll We have also seen that the 
angular momentum J tends to be uniform in the interior due to the poloidal circulation, as is clearly shown by the 
central void in the contour plot of J in Fig. I12f d) . 
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FIG. 9: Comparison of v v profiles obtained by simulation and experiment. 
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FIG. 10: Time variation of v v profile for Re = 1600. The rigid rotation profile of the endcaps and outer cylinder shown for 
comparison. 



The structures of the azimuthal and poloidal flows are summarized schematically in Fig. IT^l 

One interesting but unexpected finding of these simulations is the existence of the jet. Its characteristics are 
summarized as follows: (1) The jet becomes thinner with increasing Re, and its width is similar to that of the 
boundary layers on the caps; and (2) the jet is steady and symmetric at low Re but flaps vertically above a critical 
Reynolds number between 400 and 800. 

Experimentally, it is not straightforward to confirm the existence of a jet- like radial flow directly. Since the maximum 
radial velocity in the simulation is only a few percent of the large azimuthal flow, the streaks on the camera images 
would rotate only a few degrees, too little to be resolved by our measurements. At lower rotation speeds, we were 
able visually to follow relatively large and neutrally buoyant particles in the water. We observed rather rapid outward 
motions at the midplane after the particles were "sucked" into the boundary layers at the top or bottom endcaps. 
This is consistent with the jet-like flows indicated by the simulations. 
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FIG. 11: Time variation of angular momentum J = rv v for Re — 1600. 



F. Boundary Layers 

The width of the boundary layers (including the jet) depends on Re. As Fig. H3 indicates, higher Re causes thinner 
boundary layers. Actually, the inward (negative v r ) flow on the top and bottom boundaries is induced by the same 
mechanism as classical Ekman circulation: a viscous reduction of v v in the boundary layer, leading to an imbalance 
between outward centrifugal force and the pressure gradient. The boundary layers at the endcaps in our system do 
not have uniform width. As Figs. ll2f a1. (b), and (d) indicate, these layers are thick near the inner cylinder, reaching 
roughly 10% of the vessel's height (H). The width monotonically declines with increasing r, disappearing into the 
rigid rotation part of the outer part of the fluid. [See also Fig. H5T bl .] 

For small departures from rigid rotation, the Ekman layer thickness is Se = wnere v i s the kinematic 

viscosity. Our system is very far from rigid rotation, so it is not immediately clear what to substitute for fi. If one 
uses the mean frequency then for our geometry, Se/H « 1.24i?e -1 / 2 , hence s=s 3% at Re = 1600. In fact, we 

estimate from our simulations that the fractional thickness of the boundary layers is ~ 10% at this Reynolds number. 
For small departures from a differentially rotating state, however, we believe that it is more appropriate to scale Se 
with the epicyclic frequency, 



1 dJ^ 1/2 



Or 



(19) 



This is the maximum frequency of small axisymmetric motions (inertial oscillations) in the inviscid interior of the 
fluid, so it represents the inertial forces that must be balanced by viscous ones in order to drive a radial flow along 
the boundaries. Since k reduces to 2f2 for rigid rotation, we take Se = \f2v/R. A characteristic value for n is 



4^2 \ 1/2 

* = 2( ' 2 "rT 1 ) (20) 



This leads to Se/H — 3.39i?e -1 / 2 , or ~ 8.5% at Re = 1600, which is about three times larger than the previous 
estimate and closer to the results of the simulations. 

A prediction of the latter scaling is that the Ekman-layer thickness should increase along a sequence in which 
(rffii) / (t\£1{) approaches unity (from above) while the mean rotation yTJaOi is constant. 



G. Flow in Shorter Cylinder 



The numerical simulations presented so far have been performed for cylinders of height H = 10cm, hence about the 
same as the width of the gap (r2 — r\ = 11.1cm), as in our laboratory experiment. In order to elucidate the effects of 
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FIG. 12: Flow under enforced symmetry at Re — 1600: (a) poloidal flow; (b) contours of azimuthal velocity v v \ (c) angular 
velocity fi; and (d) angular velocity J. 



the top and bottom endcaps on the fluid motion, we have also performed numerical experiments for shorter heights: 
H = 8,6,2, and 1cm. Figure ITU shows the corresponding flow profiles after nonlinear saturation. The Reynolds 
number, which is based on the cylinder radii rather than H [eq.QJ] is 1600 in all cases. Compare these with Fig.JSfe). 

In the case of the shortest height (H = 1cm), most of the fluid rotates at the angular velocity of the endcaps. This 
is not surprising, since in the limit of infinitesimal height, the fluid would "adhere" to the endcaps. The more rapidly 
rotating inner cylinder influences the flow over a radial distance comparable to H. Comparing the panels of Fig. 1141 
we see that the domain of rigid rotation shrinks as H grows, but it still exists when H = 10cm, as summarized in the 
schematic Fig. lCT b). The existence of such a region in a short Couette flow has also been reported by Dunst^. 

H. Spin Down 

To better quantify the effects of circulation, a series of experiments and simulations have been performed to study 
the transient flow when both cylinders (and endcaps) are suddenly stopped. The rate of spin down reflects the 
efficiency with which the circulation transports angular momentum and the viscous coupling to the walls. In the 
experiments, starting from steady rotation both cylinders were braked to a complete stop within about one second. 
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FIG. 13: Schematic flow structure, (a) The Poloidal flow in the boundary layers and the jet. (b) Azimuthal flow structure. 



The flow speed in a small volume was measured against time, as shown in Fig. 1151 

A simple exponential fit to the measured data is not appropriate because the spin-down time r itself depends on 
angular velocity: 

H " (21) 



where the factor 2 comes from the fact that the circulation has 2 cells and the Ekman layer thickness Se is taken to 
be yj v/Vt (f2 is an averaged angular velocity). Thus, we have 

^ = _^ oc _03/2 ; (22) 

at t 

which leads to 

m - j^—,- (23) 

The measurements are fitted to ea. H23[) where the steady state angular velocity O(io) is known while to and r are 
fitting parameters. The fitted line is shown by the dashed line in Fig. 1151 and the spin down time, t = 11.2 ± 0.9 sec, 
is obtained. 

Spin down has also been simulated by imposing a sudden stop of all boundaries. Figure ^] shows results for 
azimuthal velocity in the same volume as in the experiment for Re — 3200. Again, the spin down time is obtained 
by fitting ea. l|23|) . with the result r = 0.82 sec. Figure IT7I displays r determined similarly for a series of simulations 
at different Re. The trend is well fit by a power law (dashed line), 

t = 0.012i?e 053 sec. (24) 
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For the purposes of Fig. 1171 and ea. (|24[l . we fix Oj = 2000 rpm and = 150 rpm (as in the experiment) and imagine 
that the Reynolds number of the simulations is controlled by varying the viscosity. If the viscosity is fixed and the 
rotation rates vary (as would be more convenient in an experiment) then the spindown time scales as Re~ 0AJ . The 
power law (PU) agrees excellently with the simple estimate given by ea. lf^Tl) . 

t = —= = — - = 0.011i?e 1/2 sec. (25) 

2VM 2y/r 1 (r 2 - n^flRe- 1 

[For comparison, ea. (|A9|) predicts = 0.0095-Re 1 / 2 sec] The experimental point, which was not included in the fit, is 
rather close to the extrapolation of ea. (|24|l : 11.2 ± 0.9 sec observed vs. 17 ± 0.9 sec predicted, using v = 0.01 cm 2 sec -1 
for water. 
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FIG. 15: Experimentally measured azimuthal velocity in the region (r, z) £ (11.5 ± 1, 3 ± 0.5) cm during spin-down after both 
cylinders and endcaps are stopped at approximately t = — 1 sec. Dotted line is azimuthal velocity in steady state when Sli = 
2000 rpm and ^2= 150 rpm. Dashed line is a fit of the form 12311 . 
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FIG. 16: Like Fig. 1151 but for a simulation at initial Re — 3200. The boundaries are stopped at t = 6.0747 sec. 
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FIG. 17: Spin down time versus Re for simulations (asterisks) and experiment (diamond). Dotted line is a fitted curve to the 
results from simulation only. 



V. SUMMARY AND DISCUSSION 

Motivated by our proposed magnetohydrodynamic experiment 0, we have studied a short, wide-gap, hydrody- 
namic Couette flow by experiment and numerical simulation. A challenge of the gallium experiment will be to set up 
an initial rotation profile that is stable to the Taylor-Couette instability (TCI) while unstable to the magnetorota- 
tional instability (MRI) when an appropriate magnetic field is applied. It is found experimentally that the azimuthal 
flow Vtplr) is significantly different from that in an infinitely long circular Couette system. In order to understand 
the underlying physics, numerical simulations have been performed using the same boundary conditions as in the 
experiment. The measured profile of v v is successfully reproduced by the simulations, which show a strong poloidal 
circulation driven by boundary layers at the endcaps. Furthermore, excellent agreement between experiment and 
simulation is found for the spin-down time when the Reynolds number is scaled. 

These agreements are rather remarkable considering that there is a factor of 300 between the Reynolds numbers 
of the simulations and of the experiments, and that the simulations are performed in two dimensions, while the 
experiments are three-dimensional. The suggestion is that the essential dynamics of the system do not change when 
Re is raised from 3200 to about 10 6 . 

Is it plausible that the the boundary layer remains laminar up to the experimental Reynolds number Re = 8.8 x 10 5 ? 
Perhaps: nonrotating flow over a flat plate remains laminar below -Re cr ;t ~ 5 x 10 5 ; perhaps more relevantly for our 
experiment, -Re cr ;t ~ 3 x 10 5 for a disk spinning freely in an extended, nonrotating fluid, where Re = SIR 2 /v based on 
the angular velocity and radius of the disk[23. For fully turbulent boundary layers, the stress (a — lateral force per 
unit area) exerted on the boundary is parametrized by a friction coefficient Cf = a/ \pV^, where Vao is the relative 
velocity of the fluid well outside the boundary layer. It is known that Cf varies slowly with Re when Re is large. 
In flow over smooth surfaces at Re ^S> 10 6 , for example, von Karman prescribes 29] Cf — 0.455(log 10 Re)~ 25S . Let 
us suppose that we can take Cf to be constant over the entire boundary (cylinders and endcaps) for the purpose of 
estimating the spin-down rate. Taking advantage of the fact that the specific angular momentum varies slowly within 
our steady-state Couette flow and defining J = (rf Cll+r^Sl 2 ,)/^, we estimate that the total torque on the fluid shortly 
after the cylinders stop is T w 2np{H + r 2 — r\)J 2 Cf. The total angular momentum of the fluid is L w ^p{r\ — r\)HJ . 
Thus the spin-down time becomes 

Tf = i = 7^ ^TTW^ 7 ~ 0.015C7 1 sec. 

1 T {r 2 n 2 + r 2 tt 1 ){H + r 2 - n) s 

This agrees with the observed value of Usee for Cf = 1.3 x 10 -3 . For comparison, von Karman's formula predicts 
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C/(10 6 ) ps 4.5 x 10 3 . Our Reynolds number however, is not far from i?e cl it, so turbulence may not be fully 
developed. Indeed, alternative definitions of the Reynolds number fall even closer to the critical value: for example, 
J/v re 3.0 x 10 5 . 

Detailed analysis of the simulations show that the poloidal circulation consists of two cells. A strong radially 
inward flow forms near each of the endcaps in a thin boundary layer. After turning into a vertical flow along the inner 
cylindrical wall, these layers merge at the midplane into a jet-like, radially outward flow to complete the circulation. 
The existence of such a jet-like feature appears not to have been recognized previously. Dunst performed a water 
experiment in short cylindrical annulus with similar condition to ours[22j. One set of five experiments by Dunst was 
carried out with rigid endcaps fixed to the outer cylinder. (Other experiments were done with a free upper surface.) 
Dunst reported the formation of a two-cell pattern as well as a region of rigid-body rotation in the outer part of the 
fluid. However, there was no description of a jet between the cells. It is possible that the jet was just overlooked. 
We note that although we have preliminary evidence for a jet-like flow at the midplane between two cells in the 
experiment, it is difficult to visualize it and to measure its detailed characteristics. We also carried out simulations 
in which the outer cylinder is stationary, leading presumably to TCI instabilities. It is interesting that the jet does 
not form in the latter situation, which has been the main focus of experimental effort in short Couette flows. 

We have seen that a region of rigid-body rotation occurs in the outer part of the system. This has been explained 
by the tendency of the flow to "adhere" to the outer cylinder and endcaps. An alternative setup of the apparatus 
would be to have the caps rotate rigidly with the inner rather than the outer cylinder. In that case, rigid-body 
rotation would be expected to appear in the inner part of the flow. Actually, such a fluid-dynamical system has been 
investigated in the literature in connection with hard disk drives for computers [3(1 l3ll l32l l33l l34| . It has been shown 
in these studies that most of the fluid rotates rigidly with the disks. 

The poloidal circulation, and especially the jet, found in this study are an interesting phenomenon in rotating 
fluids. They transport angular momentum efficiently and reduce the free energy available for shear-driven instabilities. 
Therefore, we will need to minimize this circulation in the MRI experiment. 

One idea for reducing the effect of endcaps is to use a tapered section[35|. Another idea is to divide the cap into 
two parts: the inner one fixed to the inner cylinder, and the outer one fixed to the outer cylinder. Here, we expand 
on the latter idea. We divide the endcaps into multiple rings that rotate independently (see Fig. 1181) . The angular 
velocity of each ring is chosen according to its center radius and the circular Couette flow. 

Figure shows the profiles of v v at Re = 1600 for three choices of the number of rings. The parallel contour lines 
in the figure indicate that the fluid is in a Taylor-Proudman state. The effect of the endcaps is highly localized in 
its vicinity and the poloidal circulation is suppressed almost perfectly. Without the circulation, the stationary flow 
in this differentially rotating system is very close to that of an infinitely long, ideal circular Couette flow. This is 
confirmed by the profile at the midplane when 3 and 5 rings are used (Fig. 1201) . 

The rings will reduce the relative velocity between the boundary and the interior flow in steady state, to the point 
where a turbulent boundary is unlikely (see the discussion above). This makes it more likely that the simulations 
accurately predict the interior flow despite their relatively low Reynolds number. 

An optimization process incorporating this idea is underway to design a gallium experiment with maximum con- 
trollability of the flow profile, yet with minimum complications to engineering and experimental operations. 
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FIG. 19: Simulated azimuthal velocity v v for Re = 1600 when the endcaps are divided into multiple rings; (a) 3 rings; (b) 4 
rings; and (c) 5 rings. 
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By similar steps, ea. <|13|) yields 



APPENDIX A: EKMAN LAYER IN A DIFFERENTIALLY ROTATING FLOW 

Consider a steady axisymmetric flow that departs only slightly from a centrifugally stable Couette profile, 

fl = a + -j, a, b > & constant. (Al) 

One may linearize the equations of motion 112() - (|16|l about the basic state (|A1|) . Marking first-order quantities with 
the prefix 5, we have ip — > Sip, 0Jtp — > Su> v , and v v — > rCl + Sv v . Linearizing ea. p2|l and discarding time derivatives, 

2Qd z dv v + i/Alu v = 0, 

where A = (V 2 — r~ 2 ). Incompressibility implies d z Sv z = —r~ 1 d r (r8v r ), whence 

d z 5uj v = d z (d z Sv r — d r 5v z ) = ASv r , 

so that 

2nd 2 z Sv v = -vA 2 8v r . (A2) 
1 d J r . 

— — ov r = vl\v v . (Ao) 
r dr 

For the Couette profile IjAlj) . the coefficient of Sv r in this last equation is 2a, a constant. Eliminating 5v r between 
eas.l|A2 )l &;l|A3 )l yields 

( K 2 d 2 z +v 2 A 3 )Sv v = 0. (A4) 

When viscosity can be neglected and k 2 > 0, eas. (|A2fl - (|A4|l imply that small, low-frequency n) motions are 
independent of z: this is the Proudman theorem. 

We apply these equations to the boundary layer at one of the endcaps, which rotates differentially (as might be 
achieved by dividing it into many rings, see Jvj) with a slightly different profile fl ca , P (r) from that of the fluid in the 
interior, ea. l|Alf) . Vertical derivatives are much larger than radial ones in the boundary layer, so ea. !A4fl implies that 
the boundary-layer behavior is Sv^p cx exp(A:^) with 

k A = -U/v) 2 , k = (±l±i)*P^. (A5) 

V 2v 

Of these four roots, only the two for which Sv v decays with distance from the boundary are admissible. To be definite, 
let us consider the lower endcap, so that 5R(fc) < 0. The Ekman layer thickness is 



21 



and 



Sv v = r(n cap -n)e- z ^ E cos(z/5 E ) 7 



'-cap 

Sv r = -(n cap ~n)e- z/SE sm(z/d E ) (A7) 
a 

if we take z — at the endcap rather than the midline of the cylinders. The radial mass flow is 

oo 

M E = 2irpr J Sv r dz = ■npr 2 (Q cap - SI) - 8 E . (A8) 
o 

The net torque exerted on the fluid by both endcaps is 

T = 2 J M E ^-dr = Anp J kS e r 3 (O cap - CI) dr , 
n ri 

and the amount of angular momentum that must be added to the fluid to make its rotation profile agree with the 
endcaps is 

L' — L = 2npH J r 3 (Sl cap -Q)dr. 
ri 

We may estimate the spin-up or spin-down time as 

L'-L H , A „. 
t E = « -=• (A9) 

Here k is a weighted average over radius; if we approximate it by ea. (|20|l then ea. (|A9|) predicts t E ^ 9 sec for water 
with the values of r%, ri, Cli, CI2, and H in Figs. (JIJ & l]15p. The agreement with the measured value (11.2 ± 0.9 sec) 
is perhaps better than we deserve in view of the crudeness of the theoretical treatment. In particular, since Jl C ap = 
in the spindown experiment, our linear approximation is not applicable. 
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